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SUMMARY 


Input-output characterization of fiber composites is 
studied theoretically by tracing SH waves in the media. A 
fiberglass epoxy composite is modeled as a homogeneous 
transversely isotropic continuum plate. 

The reflection of a SH wave at a stress-free plane 
boundary in a semi-infinite transversely isotropic medium is 
considered first. It is found that an incident SH wave 

reflects only a similar SH wave back into the medium. It is 
also established that the angle of reflection of the 
reflected wave is equal to the angle of incidence of the 

incident wave. 

The phase velocity of the SH waves and the delay time of 
the SH waves in reaching the receiving transducer are 

computed as functions of a reflection index, defined as the 
number of reflections of the SH waves from the bottom face of 
the continuum plate. The directivity function corresponding to 
the shear stress associated with the SH waves in the 

continuum plate is also derived as a function of the 
reflection index. 

A theoretical output voltage from the receiving 

transducer is calculated for a tone burst (a periodic input 
voltage of finite duration) . The output voltage is shown for 
tone bursts of duration 60 microseconds and center 
frequencies 0.75, 1.0 0 and 1.2 5 MHz. The study enhances the 

quantitative and qualitative understanding of the 

nondestructive evaluation (NDE) of fiber composites which can 
be modeled as transversely isotropic media. 
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INTRODUCTION 


Certain fiber composites are often modeled as equivalent homogeneous solids [1]. 
Continuum models of this sort are useful in analyzing wave propagation in fiber 
composites when the wave lengths under consideration are long compared with the mean 
fiber diameter. 

The input-output characterization of a homogeneous transversely isotropic elastic 
plate is investigated by tracing SH waves. Following the work in [2], the reflection of an 
SH wave at a stress-free plane boundary in a semi-infinite transversely isotropic medium 
is considered first. It is reestablished that an incident SH wave reflects only a similar SH 
wave. It is also reestablished that the angle of reflection is equal to the angle of incidence 
whenever the plane boundary where the reflection occurs is parallel to the isotropic plane 
of the transversely isotropic medium. 

The plane of isotropy of the equivalent transversely isotropic contiuum plate lies in the 
midplane of the plate and is parallel to each face of the plate [1]. The SH waves 
experience multiple reflections at each face of the plate. The phase velocity of the SH 
waves and the delay time for the SH waves arriving at the receiving transducer are 
computed as functions of a reflection index, which is defined as the number of reflections 
from the bottom face of the continuum plate. The directivity function corresponding to 
the shear stresses associated with the SH waves in the continuum plate is derived also in 
terms of the reflection index. 

Finally, a theoretical ouput voltage from the receiving transducer is calculated for a 
tone burst (a periodic input voltage of finite duration). The theoretical output voltage is 
expressed as an infinite sum over the reflection index. 
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REFLECTION OF INCIDENT SH WAVE AT STRESS-FREE PLANE 


BOUNDARY IN SEMI-INFINITE TRANSVERSELY ISOTROPIC MEDIUM 
WITH PLANE OF ISOTROPY PARALLEL TO PLANE BOUNDARY 

1. REFLECTED SH WAVE 

A plane progressive SH wave may be represented as 

(u,v,w)=A (P x , P y , P z ) exp{i' co(S^r + S y y +S z z -t)} (1) 

for example, see [3], where u, v and w are the displacement components of a point in the 
medium along the x, y and z axes, respectively; A is the amplitude of particle 
displacement; P x , P y and P t are the components of the unit vector of particle displacement 
along the x, y and z axes, respectively; i -V^T; © denotes radian frequency; S x , S y and S t 
are the components of the slowness vector, which points in the same direction as the 
normal to the wave front and has magnitude equal to the reciprocal of the magnitude of 
the phase velocity, see [4], along the x, y and z axes, respectively; and t denotes time. 

A plane progressive SH wave is incident on the plane boundary of a semi-infinite 
linearly elastic transversely isotropic continuum with plane of isotropy parallel to the 
plane boundary. Define a cartesian coordinate system (x, y, z) as follows: the plane 
boundary of the medium contains the x and y axes, and the z axis is the zonal axis of the 
medium (see Fig. 1). The generalized Hooke’s Law, when written relative to the (x, y, z) 
coordinate system is [4] 
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= C n ujc + C l2 v,y + C 13 w,z 
Xyy — Ci 2 u>x + C n v ,y + CjjW ,z 
x„ = Cj 3 m ,x + C n v ,_y + C 33 W ,z 

=C, 44(«^ + W ^) 

t, i = C 44 (v,z + wj) 

t^ y =C 66 («,y+v^) (2) 

where for / = j x i} is a normal stress and for i * j x tj is a shear stress; denotes partial 
differentiation with respect to the succeeding variable; C„, C I2 , C 13 , C 33 , and C u are the 
five independent elastic constants for a linearly elastic transversely isotropic medium; 
and C m is equal to 1/2 (C n -C 12 ). 

The stresses associated with a plane progressive SH wave can then be evaluated by 
substituting the expression for the displacement components at a point in the medium, 
given by eqn. (1), into eqn. (2), as 

X a = /oi4(CjjiS' jr P x + C l2 S y P y + C l3 S z P z ) exp{i o)(5p; + S y y + S 2 z — r)} 

Xyy = i(nA(C n S x P x + C n S y P y + C l3 S z P z )exp{ia(S x x + S y y+S z z-t )} 
x„ = j coA (C l3 S x P x + C l3 S y P y + C 33 S Z P X ) exp{/ co(^ + S y y + S z z -t)} 
x a = i(&A ( C 44 S t P x + C^SJ 3 ,) exp{i £0(5^ + S y y + S,z - r)} 

Xy, = i(&\(C M SJP y + C iA S y P z ) exp { i 03(5^ +S y y + S,z - r)} 

x yx = i coA (C^S y P x + C^JPy) exp{i 0 )(Sp: + S y y + S z z - 1 )} (3) 

The stress boundary conditions on a stress-free plane boundary require [3] 
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*^ = 0 

< ) +< ) =0 

^+tf = 0 (4) 

where the y> represent the stresses associated with the incident wave and the 
represent the stresses associated with the reflected waves. The stresses determined by 
eqn. (3) satisfy the boundary conditions, eqn. (4), only when the frequency, cd, of any 
reflected wave is equal to the frequency of the incident wave [3] and 

s( n = s m 

Sf = S { * ) (5) 

where Sf> and S y <’> are the x and y components of the slowness vector of the incident 
wave, respectively, and Sj*> and S y <*> are the x and y components of the slowness vector of 
any reflected wave, respectively [2]. Eqn. (5) indicates that the incident and the reflected 
waves lie in the same plane, called the plane of incidence. Without loss of generality, the 
plane of incidence can be assumed to coincide with the y-z plane; the analysis is thus 
simplified. Then, the x components of the slowness vectors of the incident and reflected 
waves can be allowed to vanish, thereby satisfying the relation 

= sf ) = 0 (6) 

An SH wave with slowness surface in a plane containing the zonal axis of a 
transversely isotropic medium possesses a transverse displacement [5]. Thus, for an SH 
wave traveling in the y-z plane in the (x, y, z) coordinate system, the unit vector of 
particle displacement is given by 

{P x ,P y , />,) = (!, 0,0) (7) 
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A wave with P x , P y and P t given by eqn. (7) and x component of slowness vector given by 
eqn. (6) has stresses determined from eqn. (3) satisfying 

^ = ^yy=^z=Sz =0 

T **0 ( 8 ) 

A P wave with slowness surface in a plane containing the zonal axis of a transversely 
isotropic medium possesses a vanishing transverse displacement [5]. Therefore, a P wave 
traveling in the plane x=0 in the (x, y, z) coordinate system has a unit vector of particle 
displacement 

(P x ,P y ,P t ) = (0,P?\P? } ) (9) 

where Pp and Pp are, respectively, the directionally dependent y and z components of 
the unit vector of particle displacement of a P wave. 

A wave P x , P y and P t satisfying eqn. (9) and S x determined by eqn. (6) has stresses 
given by eqn. (3) satisfying 

^ = ^ = 0 

V T «* 0 (10) 

An SV wave with slowness surface in a plane containing the zonal axis of a 
transversely isotropic medium also has a vanishing transverse displacement [5]. The unit 
vector of particle displacement of an SV wave traveling in the y-z plane in the (x, y, z) 
coordinate system is then 

(P x ,P y ,P t ) = (0,Pl SV) ,Pf V} ) (11) 
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where Pj SV) and Pw are, respectively, the directionally dependent y and z components of 
the unit vector of particle displacement of an SV wave. 

A wave with unit vector of particle displacement determined by eqn. (11) and x 
component of slowness vector given eqn. (6) has stresses evaluated from eqn. (3) 
satisfying 

( 12 ) 

An SH wave is incident on the plane boundary of a semi-infinite transversely isotropic 
medium. The stresses at the point of incidence associated with the incident SH wave are 
constrained by eqn. (8) and satisfy 

x (f) = x (/) = 0 

yz v zz w 

(13) 

where the superscript m is used to denote properties associated with the incident wave. 

The stress boundary conditions for any wave reflected by an incident SH wave can be 
expressed by using eqn. (4) and eqn. (8) as 

^'=^=0 

■fi 1 — l£i» 0 (14) 

where the superscript <*> is used to denote properties associated with the reflected wave. 
Eqn. (10) indicates that the and x„ components of the stress tensor associated with a P 
wave are not zero and eqn. (12) indicates that the -q. and t„ components of the stress 
tensor associated with an SV wave are also nonzero. Therefore, no P wave or SV wave 
will be reflected back into the medium by an incident SH wave, because a reflected P 
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wave or SV wave would have associated stresses that would violate the stress boundary 
conditions. Eqn. (8) indicates that a reflected SH wave can satisfy the boundary 
conditions. Therefore, an incident SH wave will cause only an SH wave to be reflected 
back into the medium, as shown in [2]. 

2. SLOWNESS SURFACE OF SH WAVE 

The equations of motion relative to the (x, y, z) coordinate system are [4] 

X xx jc+X xy ,y+T„,z=pu,tt 
x xy pc+x yy ,y + x y ,,z = pv,tt 

+x yi ,y + x„,z = p w,tt (15) 

where p is the density of the medium and ",tt" denotes the second derivative of the 
proceeding variable with respect to time. 

When the components of stress are calculated according to eqn. (3) and the 
displacement components are calculated according to eqn. (1), the equations of motion 
can be written as [2] 

(C, A 2 + C(J>y + C^S 2 - p )P X + (C 12 + CJSJSyP, + (C 13 + CJSJf, = 0 
(C 12 + CJS x S y P x + (C^l + C n S 2 y + C^S 2 - p)P y + (C 13 + CJS y S z P, = 0 
(C 13 + CJSJSf, + (C 13 + CJS y S t P y + (C M S 2 + C„S 2 y + C 33 S 2 - p)P, = 0 (16) 

Eqn. (16) can be rewritten in matrix form as 

[R](^^) r = [0] (17) 

where [fl] is a 3x3 matrix with entries 
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b n = (C n S x + C(Js y + C 44 S I — p) 

= (^66^x + Cll^y ^44 ~ P) 

633 = (C 44 S X + C^Sy + C 33 S X — p) 
bn ~ b 2 \ — (Q2 + C(^)S X S y 
b\3 = ^31 = (Cl3 + C 44 )S x S l 

b^h = ^32 = (^13 "b C M )S y S z ( 18 ) 

(PJPyf) is the unit vector of particle displacement and [0] is the 3x1 zero matrix. The 
matrix [B] as defined by eqn. (18) is a symmetric matrix; therefore, there exist three real 
eigenvalues. One eigenvalue corresponds to an SH wave, another eigenvalue corresponds 
to a P wave, and the third eigenvalue corresponds to an SV wave [2]. 

The plane wave solution is found by setting the determinant of matrix [£], the matrix 
of coefficients of P x , P y and P t , equal to zero [6]. Expanding the determinant of [5] and 
solving for the three roots yields equations for the three slowness surfaces. The slowness 
surface for an SH wave is [7] 

Cntf+S^ + C^S^ p (19) 

The values of the material constants C m , C u and p are given for a representative 
fiberglass epoxy composite in [4] as ^=3. 243x10 s N/m 2 , C w =4.422xl0 9 N/m 2 and p 
=1850 kg/m 3 . The slowness surface of an SH wave traveling in the y-z plane in the 
fiberglass epoxy composite is obtained by substituting the numerical values of the 
constants into eqn. (19) and setting the x component of the slowness vector equal to zero. 
The first quadrant of the slowness surface thus obtained is shown in Fig. 2. 
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3. ANGLE OF REFLECTION OF REFLECTED SH WAVE 


According to eqn. (5) the y component, S y <‘> , of the slowness vector of the incident SH 
wave is equal to the y component, S y <*> , of the slowness vector of the reflected SH wave. 
The z component, Sj R > , of the slowness vector of the reflected SH wave is found from the 
equation for the slowness surface of an SH wave, eqn. (19), by setting 5 w equal to S y "> 
and Sj K > equal to zero, according to eqn. (6). Then, Sj*> is found to satisfy 

Sf ) = ±S I (/) (20) 

where S<‘> is the z component of the slowness vector of the incident SH wave. The 
slowness vector of the incident wave points out of the medium, and the slowness vector 
of the reflected wave points into the medium; therefore, the relationship between S<*> and 
S< 1 > must be [2] 

Sf ) = -S 1 </) (21) 

The angle of incidence is defined as the angle between the slowness vector of the 
incident wave and the normal to the boundary at the point of incidence. The angle of 
incidence Q m of the incident SH wave is then 


0 (/) = arctan 


fs^ 

-s? 


( 22 ) 


where the point of incidence is taken to be on the plane boundary of the semi-infinite 
body (see Fig. 1). 


The angle of reflection is defined, in a manner similar to the angle of incidence, as the 
angle between the slowness vector of the reflected wave and the normal to the boundary 
at the point of reflection. The angle of reflection 0w of the reflected SH wave is then 
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Q {R) = arc tan 


(SjV) 

v 5f> y 


(23) 


Evaluating 0w using eqn. (5) and eqn. (21) and comparing the resulting expression to 
eqn. (22) establishes 


Q {R) = 6 (/) (24) 

Thus, the angle of reflection of the reflected SH wave is equal to the angle of incidence of 
the incident SH wave [2], 
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INPlJT-OIJTPlJT CHARACTERIZATION OF FIBER COMPOSITE 

1. PHASE VELOCITY AND DELAY TIME 

Certain fiber composites may be modeled as a linearly elastic transversely isotropic 
homogeneous continua [1]. A fiber composite modeled as such a solid in the form of an 
infinite plate where the plane of isotropy lies in the midplane of the plate is to be studied. 
A cartesian coordinate system (x, y, z) is chosen such that the x-y plane coincides with 
the plane of isotropy and the plate is bounded above by the plane z = h / 2 and bounded 
below by the plane z = - hi 2. 

A transmitting transducer and a receiving transducer are assumed to be coupled to the 
top face of the plate in the y-z plane and separated by a distance L. The input electrical 
voltage of the transmitting transducer is a known function of time, V.(r). The output 
electrical voltage of the receiving transducer is an unknown function of time, V o (t). The 
transmitting transducer converts the input electrical voltage into a stress that travels 
through the plate as stress waves. The receiving transducer converts the stress associated 
with the stress waves into an output voltage. In the following analysis, only those stress 
disturbances at the receiving transducer associated with SH waves generated by the 
transmitting transducer will be considered. 

The SH waves produced by the transmitting transducer experience multiple reflections 
at each face of the plate before reaching the receiving transducer. Since the isotropic 
plane of the plate lies in the midplane of the plate and is parallel to the top and bottom 
faces of the plate, each reflection may be treated as the reflection of a plane progressive 
SH wave on the plane boundary of a semi-infinite linearly elastic tranversely isotropic 
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continuum with its isotropic plane parallel to the boundary [2]. Therefore, at each 
reflection the incident SH wave produces only a reflected SH wave and the angle of 
reflection of the reflected SH wave is equal to the angle of incidence of the incident SH 
wave. 

The transmitting transducer and the receiving transducer are coupled to the same face 
of the plate, therefore, any SH wave traveling from the transmitting transducer to the 
receiving transducer must experience an odd number of reflections. Let the total number 
of reflections experienced by an arbitrary SH wave traveling from the transmitting 
transducer to the receiving transducer be equal to 2n-l, where n is called the reflection 
index and is a positive integer equal to the number of reflections from the bottom face of 
the plate. Assume the angle of incidence of the SH wave at the receiving transducer is 9,. 
The angle of reflection is equal to the angle of incidence at each reflection experienced 
by the SH wave, so the angle of incidence at each reflection is 0„ (see Fig. 3). The 
distance r a traveled by the SH wave between each reflection is calculated from Fig. 3 as 


h 

COS0„ 


(25) 


The distance b H traveled in the y direction by the SH wave between each reflection is 

b H =r n sin0 B (26) 

Using eqn. (25) to evaluate r n in eqn. (26), the distance b n can be expressed as 

= /i tan 0 n (27) 
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The SH wave travels a distance b „ in the y direction between each reflection, between 
the transmitting transducer and the initial reflection, and between the final reflection and 
the receiving transducer, so the total distance B n traveled in the y direction by the SH 
wave is 


B n = 2 nh tan 0„ (28) 

The transmitting transducer and the receiving transducer are separated by a distance L in 
the y direction, therefore; the total distance traveled by the SH wave must be L. Eqn. (28) 
can then be expressed as 

L = 2nh tan 0„ (29) 

The angle of incidence at the receiving transducer must therefore satisfy 


0„ = arctan 



(30) 


where n is a positive integer. 

Eqn. (30) allows all the trigonometric functions of 0, to be calculated (see Fig. 4). The 
tangent of 0, is given by 



(31) 


The sine of 0. is given by 



and the cosine of 0, is given by 


(32) 
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2 nh 

cos 0„ = , 

^L 2 + 4n 2 h 2 


(33) 


The distance r, traveled by the SH wave between each reflection can be evaluated 
using eqn. (25) and eqn. (33) as 


r n =‘ 


VL 2 +4n 2 /t 2 
In 


(34) 


The SH wave travels a distance r„ between each reflection, between the transmitting 
transducer and the first reflection, and between the last reflection and the receiving 
transducer; so the total distance /?„ traveled by the SH wave between the transmitting 
transducer and the receiving transducer is 

R H =^L 2 +4n 2 h 2 (35) 

The phase velocity C OT (0) of an SH wave traveling in a transversely isotropic 
continuum is [4] 

C SH (0) = {C 66 sin 2 e + C 44 cos 2 0} 1/2 (36) 

Using eqn. (32) and eqn. (33), the phase velocity can be expressed in terms of n as 


^sh(®j») 



L 2 

L 2 + 4n 2 h 2 


+ C 


44 


4 n 2 h 2 
L 2 +4n 2 h 2 


1/2 


(37) 


The time delay t, is defined as the time taken for the SH wave to reach the receiving 
transducer after being produced by the transmitting transducer. The time delay is thus 


K 

" c„(0.) 


( 38 ) 
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Eqn. (38) can be evaluated using eqn. (35) and eqn. (37) as 

t„ = { L 2 + 4n 2 h 2 } {CJ} + C u 4n 2 h 2 } ia (39) 

The time delay for a typical transducer geometry is shown in Table 1. The transducers 
are assumed to be coupled to a fiberglass epoxy composite plate having continuum model 
with material constants C 44 =4.422xl0 9 N/m 2 and C 66 =3.234xl0 9 N/m 2 and separated by a 
distance L= 10 cm. The plate is assumed to be of thickness h=5 cm. The time delay is 
shown as a function of the reflection index for values of n from 1 to 9. 

2. ASSUMPTIONS ON TRANSDUCERS 

The transmitting transducer and the receiving transducer are assumed to be transverse 
transducers that transform an electrical voltage into a uniform shear stress or a uniform 
shear stress into an electrical voltage, respectively. The following approach parallels that 
of [8]. 

If an input voltage V s of amplitude V and frequency co is applied according to 

V i = V exp {-/cor} (40) 

where i - and t denotes time, the shear stress that is introduced into the medium at 
the transducer-medium interface by the transmitting transducer is 

xjt) = F j(gj)V exp{-i(©r + <}>,)} (41) 
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where F,(co) is the frequency dependent transduction ratio for the transmitting transducer 
in transforming a voltage into a stress and is a phase angle. In eqn. (43) and eqn. (44) 
the harmonic character of the signals is expressed in complex notation, but only the real 
parts of these and subsequent equations should be considered. The amplitude T of the 
applied stress is then 

T = F 1 ((o)V (42) 

Similarly, if a stress wave producing a shear stress component x„’ of amplitude T ' and 
frequency co that impinges on the receiving transducer is defined as 

x xt '(t) = T'cxp{-iM} (43) 

then the output voltage V a from the receiving transducer is 

V 0 (t ) = F 2 (co)T'exp{-/((or + <y} (44) 

where F 2 (gj) is the frequency dependent transduction ratio for the receiving transducer in 
transforming a shear stress to a voltage, and <t> 2 is a phase angle. Thus, the amplitude V of 
the output electrical voltage is 

V' = F 2 (co)r (45) 

The characteristics of F,( <n) and F 2 ( co) are unknown except that the product F,(co)F 2 (co) is 
dimensionless. 

3. DIRECTIVITY FUNCTION 

The directivity function of the shear stress, t = , due to an SH wave traveling in a 
transversely isotropic medium is evaluated from the far-field asymptotic solutions of the 
displacement components produced by a harmonic point load imbedded in an infinite 
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body in [2]. The directivity function D< SH > is given by [2] 


D {sri) = - 


C44C< 


47t(C 66 /p)(C 44 /p) ,/2 L l 


^66 
yo ’ 


\2 f C y 

vP -j 


1/2 


(46) 


where S’ and S’ are the y and z components, respectively, of a point on the slowness 
surface of an SH wave where the normal to the slowness surface is parallel to the line 
connecting the transmitting transducer and the point at which the directivity function is to 
be evaluated. 


The shear stress, radiated by the transmitting transducer, that is incident at the 
receiving transducer with angle of incidence 0,, is assumed to be equivalent to the shear 
stress at a point M \ associated with an SH wave traveling in a semi-infinite transversely 
isotropic medium (see Fig. 5). Except for the reflection coefficients at each face of the 
plate the shear stress can be computed as if there were no bottom boundary to the plate 
[ 2 ]. 

The slope, m, of the line connecting the transmitting transducer to the point M t is 
found from Fig. 5 to be 



2 nh 
~L~ 


(47) 


Eqn. (47) establishes that the line connecting the transmitting transducer to the point M n 
lies in the y-z plane. Therefore, a point on the slowness surface, (S xt \ S ym ’, S, t ’ ), where the 
normal is parallel to the line connecting the transmitting transducer to the point M n may 
be found by setting S x ’ equal to zero. 
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The normal to the slowness surface of an SH wave is found at any point in the y-z 
plane by setting the x component of the slowness vector, S x , equal to zero and 
differentiating the equation for the slowness surface, eqn. (19), with respect to S,, the z 
component of the slowness vector. Solving the subsequent equation for the slope, 
-dS JdS t , of the normal yields 


dS y _ C^S, 
dSj Cffriy 


(48) 


Combining eqn. (47) and eqn. (48) establishes the constraint on and S„\ the point on 
the slowness surface where the normal is parallel to the line connecting the transmitting 
transducer to the point Af„, as 


_ 2nh 

CJ7.~' L 


(49) 


The point (S,/, S yx , S,/) must lie on the slowness surface of an SH wave and have 
slope satisfying eqn. (49). Therefore, S yH ‘ and S I(1 * can be found by setting S,/ equal to 
zero and then solving eqn. (19) and eqn. (49) simultaneously. Then the value of S y / is 



(50) 


and the value of S,/ is 



2nh 




L 

C44 

Cu + CfA 

[&)] 


( 51 ) 
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Substituting eqn. (50) and eqn. (51) into eqn. (46) produces an expression for Dj SH >, 
the value of the directivity function associated with an angle of incidence 0, at the 
receiving transducer. The expression for Df H > as a function of n is 


p. (M) _ nh CO 
" L 2n 



(52) 


The directivity function for an SH wave traveling in a transversely isotropic plate is 
calculated as a function of the angle of incidence at the receiving transducer by treating 
eqn. (30) and eqn. (52) as a pair of parametric equations with parameter n equal to the 
reflection index. The directivity function D< SH > is shown for an SH wave traveling in a 
fiberglass epoxy composite plate as a function of the angle of incidence at the receiving 
transducer in Figs. 6, 7 and 8. The transducer geometry is assumed to be L=10 cm and 
h = 5 cm. The composite plate has with material properties (7^=4.422x1 0 9 N/m 2 and 
(7 66 =3.234xl0 9 N/m 2 [4]. The directivity function is shown for frequencies equal to 0.75, 
1.00 and 1.25 MHz. 


4. STRESS FIELD RADIATED BY TRANSMITTING TRANSDUCER 

If there were no bottom boundary to the plate and the SH waves were propagating in 
an infinite half-space, the SH wave incident at the receiving transducer with angle of 
incidence 0„ would travel to a point M,, in a time t M , traveling a distance /?. (see Fig. 6). 
The amplitude of the hypothetical shear stress at the point M, is Tj u > and is defined as [2] 

£)(SH) 

-%-expl-aR.} (53) 
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where T is the magnitude of the point load applied by the transmitting transducer in the z 
direction, Dj SH >, the value of the directivity function for an SH wave associated with the 
point is given by eqn. (52) and a is the SH wave attenuation constant of the medium. 
The SH wave, however, does not propagate in an infinite half-space, but instead 
experiences 2n-l reflections. The amplitude T 0it of the shear stress at the receiving 
transducer is thus obtained by modifying eqn. (53) as 

Tn VH) (n (SH-SH\2A-l 

Ton = - >— expt-oKJ (54) 


where Q SII - SH > is the reflection coefficient for SH waves and is equal to -1 [2]. 

The amplitude V 0n of the output voltage from the receiving transducer is obtained from 
eqn. (45) and eqn. (54) and is 


T , -F 2 (gs)TD® H) r b , 
V, H = R ex P{“«^n} 


(55) 


Substituting eqn. (42) into eqn. (55) yields 


v , m = £ exp{-c iR n } 


(56) 


where V is the amplitude of the input voltage. Introducing a possible electrical signal 
amplification factor, K , eqn. (56) can be written as 


V„ = 


-KF^F^aWD 


(SH) 


-exp{-a /?„} 


(57) 
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5. OUTPUT VOLTAGE DUE TO TONE BURST 


Assume the input voltage V. to the transmitting transducer is a periodic function of 
time, of center frequency go and duration t r Mathematically the input voltage can be 
expressed as the sum of two periodic functions of frequency go that are 180° out of phase, 
one beginning at time /=0 and one beginning at time t=t r The input voltage is then given 
by 

V ; = V exp{/Gor} U(r) - V exp{icor} U(r - r,) (58) 

where V is the amplitude of the input voltage, /■Vm , and UQt) is the unit step function 
defined as 


fO if x <0 


U(x)H 


(59) 


[1 if x > 0 

The periodic nature of the signal in eqn. (58) and subsequent equations is expressed in 
complex notation, but only the real part of the signal is to be considered. 


The voltage output V o at the receiving transducer is found by considering the two 
terms on the right-hand side of eqn. (61) independently and superposing the voltage 
output associated with each input signal. Consider the input voltage V' given by the first 
term on the right-hand side of eqn. (58) 


V,' = V exp{/cor} U(r) (60) 

The stress waves produced by a transmitting transducer with input voltage V! will have a 
periodic character with respect to time. Therefore, the amplitude of the output voltage V„ n ' 
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associated with the SH wave experiencing 2/j-l reflections will also be periodic in time 
and the output voltage V' at the receiving transducer must be found by the superposition 
of the contributions of each reflection path. 


An SH wave experiencing 2n-l reflections between the transmitting transducer and 
the receiving transducer will be in phase at the receiving transducer with an SH wave 
traveling in a semi-infinite body to a point AT, (see Fig. 5). The x a stress associated with a 
plane progressive SH wave is given by eqn. (3) and can be calculated at point Af, as 


= T m 

a M M n 



L-S. 2nh - 1 

x n 


(61) 


where T Mh is the amplitude of the stress and S yK and S ln are the y and z components of the 
slowness vector of the SH wave, respectively. The shear stress at the receiving 
transducer is then 



L-S , 2nh 


-)]) 


(62) 


where T 0it is the amplitude of the stress and is given by eqn. (54). The output voltage V 0n ' 
at the receiving transducer is found by substituting eqn. (62) into eqn. (44) to be 


V 0n \t) = F 2 (co)r 0)i exp|t[co^ L -S t 2nh + (63) 

where 0 is a phase angle. The output voltage is then calculated using eqn. (42) and eqn. 
(54) and including a possible electrical signal amplification factor K as 
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-KF l ((a)F 2 (Gi)VDi SH) I r ( \ l] 

V = - exp{-ctf „} expj ij© ^ L - ^2*/* - 1 J + <|»J \ (64) 

where V is the amplitude of the input voltage and the D< SH > is given by eqn. (52). The total 
output voltage V' is found by summing eqn. (64) over all n and retarding each 
contribution by the delay time t n to be 


V/= I 

n = 1 


- -KF x «£>)F 2 (a)VD? H) 


R. 


exp{-cd?„}x 



S y L-S, 2nh 

x n 



m-o 


(65) 


where t H is given by eqn. (38). 

The output voltage V" associated with the input voltage characterized by the second 
term on the right-hand side of eqn. (58) is found similarly to be 


„„ - JrF,(<B)F ! (0))V'Dr l , 

K = 2 ^ exp{-o (RJ X 


n - 1 


R. 


expj/ co js >( L - S,2nh - r j + <]) j| U(r - (r ( - - rj) 


where t. is the duration of the input signal. 


( 66 ) 


In order to evaluate the infinite sum in eqn. (65), S yn and must be found as functions 
of n. The components of the slowness vector (0, S yn , S, H ) associated with the SH wave 
experiencing 2n-l reflections between the transmitting transducer and the receiving 
transducer are constrained by eqn. (22) to satisfy 
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( 67 ) 



The output voltage V g associated with the input voltage V. determined by eqn. (58) is 
found by adding eqn. (65) and eqn. (66). Substituting in the values of the constants given 
by eqns. (35), (39) and (69) allows V 0 to be written as 
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- f D (SH) I 

V 0 (t) = —KF l (ca)F 2 (<i))V I -■ : exp{-oWL 2 + 4n 2 /! 2 } x 

n = 1 ['vL 2 + 4/r/z 2 




’ 

/ 




L 2 

exp- 

/- 

£0 

2nh 




\ 


\l/2 




+ C 


44 


+ 2nh 




+ C, 


44 


Y' 2 " 

■N 


-t 

+ 4) 

► 





R- 


{L 2 + 4« 2 /i 2 } {C^ 2 + C 44 4« 2 /i 2 }" 1/2 1 - 




Ulf-I {L 2 + 4n 2 h 2 } {CJ. 2 + C«4n 2 h 2 } xn + 1, 


i 


(70) 


where DJ SH > is a function of n according to eqn. (52). 

The output voltage predicted by eqn. (70) is shown as a function of time in Figs. 9, 10 
and 11. The transducers are assumed to be coupled to a representative fiberglass epoxy 
composite plate. The composite plate can be modeled as a transversely isotropic 
continuum plate with material properties p=1850 kg/m 3 , (7^=4.422x1 0 9 N/m 2 and 
C 66 =3.234 x 10 9 N/m 2 and thickness h=5 cm and transducer separation L=10 cm [4]. The 
continuum plate is assumed to have SH wave attenuation constant a equal to zero. The 
normalized output voltage VJKFJF 2 V is shown for center frequencies equal to 0.75, 1.00 
and 1 .25 MHz. and a signal duration of 60|is; the phase angle <j) is set equal to zero. 
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CONCLUSION AND DISCUSSION 


The input-output characterization of certain fiberglass epoxy composite plates can be 
studied using a transversely isotropic continuum plate model when the wavelengths under 
consideration are long compared to the mean fiber diameter. For wavelengths close to or 
less than the mean fiber diameter, a continuum plate model is no longer appropriate and 
the inhomogeneities of the composite must be considered. 

In the analysis of a transversely isotropic continuum plate, the angle of reflection of a 
reflected SH wave is equal to the angle of incidence of the incident SH wave and only an 
SH wave is reflected by the incident SH wave for reflections at the top or bottom face of 
the plate. This is due to the fact that the isotropic plane of the plate is parallel to the two 
faces of the plate. When this geometry is present, the plate can be modeled as a 
semi-infinite transversely isotropic medium, except for the reflection coefficients 
encountered at each reflection which must be factored in separately. Treating the 
continuum plate as a semi-infinite body reduces the complexity of tracing stress waves 
through the medium. If this special geometry does not exist, the use of this simplified 
analysis is inappropriate and the plate must be analyzed in a more complicated manner. 

It is possible to compute the relationship between the input voltage to the transmitting 
transducer and the output voltage from the receiving transducer in terms of the material 
constants, the geometry of the transducer arrangement and a single variable, the reflection 
index n. Using the fact that the angle of incidence of an incident SH wave is equal to the 
angle of reflection of the reflected SH wave, the angle 9 with which the final wave of any 
series of waves from the transmitting transducer to the receiving transducer is incident at 
the receiving transducer can be computed as a function of the reflection index. The 
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nonzero components of the slowness vector of an SH wave traveling in a plane 
perpendicular to the isotropic plane of the plate, the phase velocity of the SH wave, the 
time delay in reaching the receiving transducer, the directivity function of the shear stress 
produced by the transmitting transducer and the output voltage from the receiving 
transducer can then be computed as functions of the reflection index and the geometric 
and physical properties of the plate. 

The output voltage from the receiving transducer, as predicted by eqn. (70), is shown 
in Figs. 9, 10 and 11 as functions of time for input voltages of 60 Jis duration and center 
frequencies of 0.75, 1.00 and 1.25 MHz., respectively. All three output voltages are 
observed to exhibit similar behavior, being characterized by simple harmonic functions 
experiencing discontinuities in amplitude. Checking the times at which these changes in 
amplitude occur against Table 1, which shows the time delays for reflection indices from 
1 to 9, reveals that the amplitude changes occur at the time delays for these values of the 
reflection index. Thus the discontinuities in the output voltages correspond with the 
arrival times of different wave paths. Since the time delay is independent of the 
frequency of the wave, these discontinuities occur at the same time in each output voltage 
regardless of the frequency of the associated input voltage. 

This study enhances the theoretical understanding of the nondestructive evaluation 
(NDE) of transversely isotropic media such as certain fiber composites. 
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Table 1 Time delay in reaching receiving transducer for SH 
wave in composite plate for reflection indices n 
from 1 to 9 . 


reflection index n 

time delay t^ 
(microseconds) 

1 

2.28 

2 

3.46 

3 

4.82 

4 

6.25 

5 

7.71 

6 

9.18 

7 

10.6 

8 

12.2 

9 

13.6 


L = 1 0 cm 

h = 5 cm 

C 44 = 4.422 x 10 9 N/m 2 

= 3.230 x 10 9 N/m2 

6 6 
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Fig. 1 Angle of incidence 9^ of incident SH wave 

/ n \ 

and angle of reflection 9 of reflected 
SH wave shown in y-z plane of (x,y,z) 
coordinate system. 
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Fig. 3 SH wave experiencing 2n-l reflections between 

transmitting transducer and receiving transducer. 
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Fig. 
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4 Right triangle from which trigonometric functions 

of 0 are derived, 
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x L 



Fig. 5 SH wave traveling distance R n to hypothetical 
point M n in semi-infinite body. 
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Magnitude of Directivity function 



/ cu \ 

Directivity function D' 'of shear stress as function of 
angle from the z axis for frequency 0.75 MHz. 


Fig. 6 
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Fig. 8 

Directivity function 
of angle from z axis 

/ Ctl j 

D l ' of shear stress as a 
for frequency 1.25 MHz. 
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Fig. 10 


Normalized output voltage as function of time for frequency 
1.00 MHz. 



Time (microseconds) 

Fig. 11 Normalized output voltage as function of time for frequency 
1.25 MHz. 
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